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:  Theoretical  calculations  are  presented  for  the  probe  absorption  spectra  of  atoms  having  magnetically 

i  degenerate  ground  states  which  are  subjected  to  weak  pump  fields  that  produce  spatial  polarization  gra¬ 

dients.  Apart  from  the  Raman  resonance  contribution  to  the  spectrum  which  has  been  observed  experi¬ 
mentally,  we  predict  that  additional  resonant  structures  should  appear  as  a  result  of  the  atomic  recoil  ^  mjm 
that  occurs  during  the  absorption  or  emission  of  radiation  by  the  atoms.  The  width  of  the  atomic  recoil  1  I 
resonances  is  directly  related  to  the  Doppler  width  associated  with  the  driven  transition.  For  sub-  *  ■* 

Doppler  cooled  atoms,  the  width  can  be  much  narrower  than  the  optical  pumping  rate  that  determines  :le< 
the  width  of  the  Raman  signal.  __  Aliro 


PACS  number(s):  32.80.Pj,  32.70.  -  n,  42.50.  -  p 


I.  INTRODUCTION 

The  field  of  sub-Doppler  laser  cooling  of  neutral  atoms 
has  progressed  to  such  a  level  that  it  is  now  possible  to 
obtain  atomic  vapors  at  very  high  density  (~1010  cm"'1) 
and  extremely  low  temperature  (~2  pK)  [1],  Since  the 
Doppler  and  transit  broadening  are  greatly  reduced  in 
these  systems,  such  cold  and  dense  samples  provide  ideal 
environments  for  carrying  out  many  nonlinear  spectro¬ 
scopic  experiments.  Moreover,  one  might  expect  some 
new  and  interesting  results  which  are  absent  or  nonob¬ 
servable  at  normal  temperatures.  For  example,  the  width 
of  the  velocity  distribution  of  atoms  at  room  temperature 
is  much  larger  than  the  photon  recoil  velocity  ftk  /m, 
where  k  is  the  wave  vector  of  light  and  m  is  the  atomic 
mass.  As  a  result,  effects  associated  with  changes  in 
atomic  velocity  upon  absorption  or  emission  of  a  photon 
can  generally  be  neglected  when  considering  the  spectros¬ 
copy  of  thermal  atoms.  However,  when  the  distribution 
width  itself  is  of  the  order  of  several  ftk  /m  as  in  the  case 
of  sub-Doppler  temperature  atoms,  such  recoil  effects 
may  not  be  negligible.  It  has  already  been  shown  that  the 
inclusion  of  recoil  effects  leads  to  new  resonance  struc¬ 
tures  in  the  nonlinear  spectroscopy  of  an  ensemble  of 
"two-level"  atoms  cooled  to  sub-Doppler  temperatures 
[2].  The  resonances  are  centered  at  8  =  0,  where  8  is  the 
frequency  detuning  of  a  probe  field  from  a  pump-field  fre¬ 
quency.  The  width  of  the  resonances  is  characterized  by 
the  Doppler  width  of  the  ensemble  ku  =kp0/m,  where  u 
\s  the  most  probable  atomic  speed  and  pQ  =  mu  character¬ 
izes  the  momentum  distribution  width.  For  many  atom- 
<^neld  interactions,  the  two-level  approximation  is  not  sa- 

E factory.  This  is  especially  true  in  sub-Doppler  laser 
oling  where  the  magnetic  degeneracy  of  the  ground- 
ite  levels  plays  a  central  role  in  the  cooling  process. 

In  this  paper,  we  extend  our  previous  calculations  to 
low  for  the  magnetic  degeneracy  of  the  atomic  levels, 
ijfhe  inclusion  of  magnetic-state  degeneracy  introduces  a 
$i:w  dimension  in  the  problem.  In  both  pump-probe  and 
Hour-wave  mixing  spectroscopy  involving  magnetically 
,  Qi  generate  atoms,  it  is  possible  to  observe  narrow  reso¬ 
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nances  centered  at  8  =  0,  whose  width  is  characterized  by 
the  rate  I”  at  which  the  ground-state  sublevels  are  opti¬ 
cally  pumped  by  the  fields  [3].  It  is  not  obvious  that  the 
recoil-induced  resonances,  characterized  by  a  width  ku 
assumed  to  be  much  smaller  than  P,  are  not  obliterated 
by  the  optical  pumping  resonances.  It  turns  out,  howev¬ 
er,  that  the  recoil-induced  resonances  persist  in  the  pres¬ 
ence  of  optical  pumping  and  lead  to  narrow  resonant 
structures  superimposed  on  the  line  shapes  obtained 
when  recoil  effects  are  neglected.  Such  narrow  reso¬ 
nances  may  have  been  already  observed  experimentally 

(4]. 

To  be  specific,  we  consider  the  absorption  of  a  probe 
field  having  frequency  ft +  8  in  the  presence  of  a  pump 
field  having  frequency  ft.  Generally  speaking,  the  probe 
absorption  is  proportional  to  ground-state  density-matrix 
elements.  In  addition  to  the  ground-state  population,  the 
absorption  also  depends  on  differences  in  population 
among  different  magnetic  sublevels,  as  well  as  magnetic- 
state  coherence.  Such  quantities  are  conveniently  ex¬ 
pressed  in  an  irreducible  tensor  basis  as  p[](g),  pi  -  i<g ). 
and  p5  ,  |  ,2(g),  where  p[](g)  is  proportional  to  the 
ground-state  total  population,  and  pp(g)  and  py(g) 
represent  ground-state  orientation  and  alignment,  respec¬ 
tively.  Under  the  low-field-intensity  approximation  to  be 
defined  below,  one  can  adiabatically  eliminate  the 
excited-state  quantities  and  obtain  a  set  of  equations  in¬ 
volving  only  ground-state  density-matrix  elements  [3.5]. 
As  a  direct  result  of  these  equations,  p[](g)  remains  con¬ 
stant  in  the  absence  of  recoil  while  p£,(g)  (K^O)  decay  at 
rates  characterized  by  the  optical  pumping  rate 
T'  =  ( T|  vl2)/[( T/2)2  + A2],  where  \  is  a  pump-field  Rabi 
frequency.  T  is  the  excited-state  population  decay  rate, 
and  A  is  the  pump  detuning  from  the  atomic  resonance 
frequency.  As  a  consequence,  one  might  expect  that  the 
nonlinear  pump-probe  signals  should  have  a  width  at 
least  of  order  T'  in  the  absence  of  recoil  or  some  other 
motional  narrowing  effects. 

To  take  into  account  recoil  effect,  one  needs  to  quan¬ 
tize  the  atomic  ccntcr-of-mass  momentum  and  write  the 
density-matrix  elements  as  p(p.p').  As  a  result  of  the 
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recoil  shifts,  the  momentum-integrated  off-diagonal  ma¬ 
trix  elements  Jp^p.p  +  Apl^p  for  Ap^O  is  not  con¬ 
served.  Assuming  the  initial  p]](p,p')  is  diagonal  in  p 
space,  then  owing  to  the  interaction  of  the  atoms  with  the 
pump  and  probe  fields,  off-diagonal  elements  of  p]](p,p'l 
are  generated,  which  are  responsible  for  the  recoil- 
induced  signals. 

The  recoil-induced  probe  absorption  or  amplification  is 
most  easily  understood  in  terms  of  Raman  transitions  be¬ 
tween  atomic  center-of-mass  momentum  states,  as  men¬ 
tioned  at  the  end  of  Ref.  [2],  We  want  to  show  here  in 
some  more  detail  that  there  exists  a  close  analogy  be¬ 
tween  Raman  process  between  internal  atomic  states  and 
Raman  process  between  external  momentum  s  ates.  As 
an  example  of  the  former,  consider  the  three-level  system 
shown  in  Fig.  1(a).  The  two  ground  states  a  and  b  are 
separated  by  some  finite-energy  splitting  6 Eha  >  0.  The 
frequencies  of  the  pump  and  probe  fields  are  given  by  il 
and  iT  =  n  +  5,  respectively.  Both  frequencies  are  as¬ 
sumed  to  be  detuned  far  from  the  resonance  frequency 
between  state  a  or  b  and  the  excited  state  c;  therefore,  the 
population  of  state  c  is  negligible,  and  states  a  and  b  are 
connected  by  stimulated  transitions  involving  both  pump 
and  probe  fields.  Assume  that  pd  >ph,  where  pa  and  ph 
are  the  populations  of  states  a  and  b.  For  6>0,  there  is 
more  absorption  of  probe  photons  and  emission  of  pump 
photons  than  the  other  way  around;  as  a  result  one  has  a 
net  probe  loss  [see  Fig.  1(a)].  By  the  same  token,  for 
5  <  0,  one  has  a  net  probe  gain. 

Now  consider  the  case  of  a  Raman-type  transition  be¬ 
tween  the  atomic  momentum  states.  Assume  that  the 
atomic  system  has  some  momentum  distribution  centered 
at  p  =  0  whose  width  is  on  the  order  of  several  fik,  as 
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FIG.  1.  (a)  Energy  diagrams  for  a  three-level  atom.  The  Ra¬ 
man  transitions  are  between  the  two  ground  states  a  )  and  h  ), 
whose  energies  are  separated  by  an  amount  bEhtl .  (b)  The 
center-of-mass  momentum  distribution  of  atoms.  Raman  tran¬ 
sitions  occur  between  momentum  states  p  and  p  +  2 fik  via  ab¬ 
sorption  of  a  probe  photon  and  emission  of  a  pump  photon. 
Based  on  conservation  of  energy,  arrow  (il  corresponds  to  b  -0 
ip  >  M  i,  and  arrow  (iil  to  b  •  0  l/i  <  fik). 


shown  in  Fig.  Kb).  The  probe  and  pump  propagate  in  the 
+  k  and  —  k  directions,  respectively.  If  the  initial  atomic 
momentum  in  the  k  direction  is  p ,  then  after  absorption 
of  one  probe  photon  and  emission  of  one  pump  photon, 
the  final-state  momentum  is p  +  2fik.  The  resonance  con¬ 
dition  of  such  a  process  obtained  from  conservation  of 
energy  involving  the  center-of-mass  degree  of  freedom  is 
given  by 

6—  — —  4m;  —0  ,  (1) 

m 

where  iok=fik1  /2m  is  the  recoil  frequency.  Assuming 
that  the  momentum  distribution  remains  unchanged  dur¬ 
ing  the  course  of  interaction,  the  resonance  width  is 
essentially  determined  by  the  interaction  time  with  the 
fields.  When  6>0,  the  initial  momentum  obtained  from 
Eq.  (1)  satisfies  p  >  —  fik:  consequently,  the  final  momen¬ 
tum  is  (p  +  Ifik  )  >  fik.  As  one  can  see,  the  initial  momen¬ 
tum  state  is  more  populated  than  the  final  momentum 
state,  resulting  in  a  net  probe  loss.  Similarly,  when  6  <0. 
the  initial  momentum  state  is  less  populated  than  the  final 
momentum  state  and  one  has  a  net  probe  gain.  When 
6  =  0.  initial  and  final  momentum  state  have  the  same 
population;  there  is  no  net  probe  gain  or  absorption.  One 
can  analyze  in  the  same  way  the  other  process  of  going 
from  p  to  p  —  Ifik  by  absorption  of  one  pump  photon  and 
emission  of  one  probe  photon  and  it  leads  to  the  same 
conclusion  as  above.  The  whole  effects  are  related  to  the 
difference  of  populations  between  two  momentum  states 
differing  by  2 fik,  therefore  one  would  expect  that  these 
effects  will  become  significant  for  the  signal  generation 
only  when  the  atomic  momentum  distribution  width  is  of 
the  order  of  several  fik.  Also  from  the  above  analysis,  it's 
easy  to  see  that  the  signal  has  a  w  idth  of  order  ku  as  one 
tunes  6  across  the  momentum  distribution.  The  magni¬ 
tude  of  the  recoil-induced  signal  is  proportional  to  the 
quantity  fik  /( pt)ku  )  —  u>k  /( ku )'.  Since  the  magnitudes  of 
the  background  signals  having  w  idth  of  order  T’  are  pro¬ 
portional  to  l/r\  the  dimensionless  parameter  that 
determines  the  amplitude  of  the  recoil-induced  signal  to 
the  resonance  having  width  f  is  of  order 

r:  **  _  {2) 

ku  p()  { ku )' 

As  will  be  discussed  below,  this  quantity  can  be  compara¬ 
ble  to  unity. 

As  a  specific  example,  we  consider  the  case  of  a 
Jg  =  1  —  2  transition  driven  by  two  pump  fields  hav¬ 

ing  a  '  and  a  polarizations  and  a  it  probe  field.  Foi 
this  type  of  pump-field  configuration,  the  atomic  distribu¬ 
tion  function  is  uniform  in  space,  owing  to  the  fact  that 
the  energy  shifts  of  the  atomic  ground-state  sublevels  are 
spaee  independent  [3].  Therefore,  any  spatial  localization 
effects  of  the  atoms  can  be  neglected.  These  effects  can  be 
important  for  some  other  situations,  such  as  with  two 
perpendicularly  polarized  pump  fields,  and  they  can  in¬ 
troduce  further  complications  into  the  calculations.  We 
shall  address  these  other  pump-field  configurations  at  the 
end  of  the  article.  This  paper  is  arranged  as  follows.  In 
See.  II,  we  derive  the  generalized  optical  Bloch  equations 
[b]  for  the  atomic  density-matrix  elements,  and  invoke 
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the  low-field-intensity  and  low-energy  approximations. 
Then,  in  Sec.  Ill,  we  calculate  the  probe  absorption  sig¬ 
nal,  which  includes  the  contribution  of  the  recoil-induced 
resonances  between  the  atomic  center-of-mass  momen¬ 
tum  states,  and  that  of  the  Raman  resonances  between 
the  atomic  ground-state  sublevels.  Finally,  in  Sec.  IV,  a 
discussion  of  our  results  and  some  other  situations  with 
different  pump-field  configurations  is  given. 


II.  GENERALIZED  OPTICAL  BLOCH  EQUATIONS 

The  atomic  level  scheme,  along  with  the  relative 
strengths  of  the  transitions  (Clebsch-Gordan  coefficients 
associated  with  the  transitions)  is  shown  in  Fig.  2(a).  The 
ground  and  excited  states  have  angular  momenta  Jg  =  1 
and  Jc—  2,  respectively.  The  incident-field  configuration, 
shown  in  Fig.  2(b),  is  as  follows:  the  two  pumps  fields 
propagate  in  the  ±z  directions,  having  a  '  and  a  polar¬ 
izations,  respectively.  The  probe  field  is  a  polarized, 
propagating  in  a  direction  opposite  to  that  of  the  a 
pump.  The  total  field  can  be  written  as 

+  ,!J'  +  6  e  e  ""  ,s" 

+  6'e.  eikz  lir‘)+c.c.  ,  (3) 

where  f,  —  +  ( 1 /v  2 )(x±fy ),  and  x  and  y  are  unit  vec¬ 
tors  in  x  and  y  directions,  respectively.  This  pump-field 
configuration  is  known  to  lead  to  sub-Doppler  cooling 
[3], 

The  Hamiltonian  for  the  atomic  system  including  the 
center-of-mass  motion  is  given  by 

2  2 

H  =  -^--+fuo  5"  \em)(em\  +  V  ,  (4) 

7  ^ 


FIG.  2.  (a)  The  atomic  level  scheme  of  a  Jg  =  1  --*7,.  =2  tran¬ 
sition,  along  with  the  Clebsch-Gordan  coefficients  associated 
with  various  transitions,  (b)  The  incident-field  configuration 
for  the  Jg  =  1  — -J,.  =  2  transition.  The  forward  and  backward 
pump  fields,  both  having  frequency  ft,  are  o  '  and  a  polarized, 
respectively.  A  probe  field  of  frequency  ft  r  5  is  a  polarized 
and  propagates  in  the  same  direction  as  the  a  '  pump. 

where  p  is  the  z  component  of  the  atomic  center-of-mass 
momentum,  and  to  is  the  atomic  resonance  frequency.  In 
the  rotating-wave  approximation,  the  interaction  Hamil¬ 
tonian  V  in  the  combined  basis  of  magnetic  substates 
| g,m)  (m  —  —  1, . .  .  ,  1 )  and  \e,m')  {m'=  —  2, ...  ,2), 
and  the  atomic  center-of-mass  momentum  state  Ip  )  can 
be  written  as 


V—y.fiX‘e  lil'  ie2,p  +ftk  )(g  l,p  |  -1 — ^=- \e\,p  +fik  MgO.p  |  4 — W  ieU.p  +fik  )  {g  —  l,p  \ 

P  v  2  v  6 

+  c  ,n,k _ 2,p -ftk  )(g  —  l,p|  +fi\'e  lll'\e  —  2,p+fik)(g—\,p\ 

r 

+  (ft\  e  liu\e  -  l,p -ftk  Xg0.pl  +fix'e  ’'|i>  -  l,p  +  fik  >  <g0,p  | ) 

+  --•(«*'  e  ,n‘\e0,p—hk  Xgl.pl  -Hix'e  '"'leO.p  +fik  Hg  l,p  i )]  +  H.c.  ,  (5) 

v  6 


where 

cr  .  er  (>' 

y  .  ~  -  _  ,  y'  =  —  _  -  , 

2v'5H  2v/5« 


p  =  .  (8) 

a.n 

(6)  where  a,/?=  \g,m  ),  k\m ' ).  Introducing  an  interaction 
representation  defined  by 


and  ri>ft  is  a  reduced  matrix  element  of  the  dipole  moment 
operator.  In  deriving  Eq.  (5),  the  property 

e  ‘  ,k:\p  )  -  Ip  t  f,k  )  (7) 

of  the  atomic  center-of-mass  position  operator  z  has  been 
used  The  atomic  density  matrix  is  expanded  as 


PnJP’P'^PaJPy)  < 

(9) 

Pltn,.,m^P<P''>=Pl!m.,m^P<P'^‘""  • 

one  can  use  the  Schrbdinger  equation  with  the  Hamil¬ 
tonian  (4)  to  obtain  a  set  of  generalized  optical  Bloch 
equations  for  plt  rf,p,p').  For  example,  one  finds 
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Pg-\,e^py)=~  Y  +  i 1 


Pg-\,e-2(P’P')-‘X*~e  g-^P'P'+fik  )] 


~iX'*e‘*'[Pe-2.e  2(P+^fc*P')~Pg-1  ,g  -!</>> p' —*k  )]  , 


where 


r= r/2 

is  the  electronic  state  coherence  decay  rate,  and 

A  =  n-w,  A'  =  n'-w  =  A  +  6  .  (11) 

We  now  invoke  a  low-field-intensity  and  a  low-energy 
approximation  [6].  In  the  low-field-intensity  limit 

|*±i«r,|A|,  (12) 

quasi-steady-state  solutions  for  the  electronic  state  coher¬ 
ence  and  excited  state  density-matrix  elements  can  be  ob¬ 
tained  in  terms  of  ground-state  density-matrix  elements. 
For  example,  the  approximate  solution  of  Eq.  (10)  is 
given  by 

Pg-l.e-2  IP,P')= - 

y+i 


Xpg_,  g_i(p,p'  +  fik ) 


Xpg_ig_l(p,p'—fik).  (13) 

When  substituting  similar  expressions  for  pge,peg,  and  pee 
into  the  evolution  equations  for  ground-state  quantities, 
one  obtains  a  set  of  equations  involving  ground-state 
density-matrix  elements  only.  Furthermore,  we  assume 
the  atomic  center-of-mass  energy  is  sufficiently  small  to 
allow  us  to  neglect  the  kinetic-energy  term 
(p2—p'2)/2mfi  in  comparison  with  |A|  (or  |A'|).  Both 
the  low-field-intensity  and  low-energy  approximations  are 
generally  valid  in  regimes  appropriate  to  sub-Doppler 
cooling,  assuming,  as  we  do  in  this  paper,  that  |  A i  >  T. 
With  these  two  approximations,  one  obtains  a  set  of 
simplified  Bloch  equations  for  the  atomic  ground-state 
elements.  They  are  given  in  the  next  section. 

Using  the  low-field-intensity  and  low-energy  approxi¬ 
mations,  we  solve  for  the  probe  absorption  coefficient  to 
all  orders  in  the  pump-field  strength  and  to  first  order  in 
the  probe-field  stiength.  As  shown  in  Appendix  A,  the 
pump-field  modified  probe  absorption  is  related  to  an 
atomic  coherence  pge ,  which  can  be  expressed  in  terms  of 
the  ground-state  density-matrix  elements  as 


where  we  have  assumed  that  y+=y_  =  y,  and  the  quan¬ 
tities  on  the  right-hand  side  of  the  above  equation  are 
defined  by 

p['  \p,p  +  2fik  )=psu(p,p  +  2fik  )e  ~'6,  ,  (15) 

where 

i 

PS<P.P')=  1  P*m.gm(P-P'>  >  <16) 

m  =  —  1 

and 

Por’  =  e  JVo'rW  +2*k  )dp  , 

pii ’=*"<6,JV.!W+2«c  )dp  ,  (17) 

Pi!  - 1 =  e  _ '8(  f  Pgig  - 1 <P,P )dp  , 
where 

PoMP')=-^-[pgi.gi(P.P')-p«-i.g-i(p.P')]  • 

j  (18) 

Pa|(P.P')=-^[Pgl,gl<P.P')+Pg  l.g-i(P.P') 

“2pg0,go(p.P'>]  • 

The  superscript  (1)  of  various  p’s  or  p  s  in  Eqs.  (15)  and 
(17)  indicates  that  these  quantities  are  proportional  to  y'* 
(in  the  resonance  or  rotating-wave  approximation,  there 
is  no  contribution  to  pge  which  is  proportional  to  y'). 
Also  in  this  paper,  unless  specified  otherwise,  a  super¬ 
script  (0)  of  a  density-matrix  element  indicates  that  this 
element  is  zeroth  order  in  probe-field  strength,  while  a 
density-matrix  element  without  any  superscript  means 
that  it  is  both  zeroth  and  first  order  in  probe  strength. 

Before  going  into  the  details  of  the  calculation,  it  is 
helpful  to  identify  various  terms  in  Eq.  (14)  with  different 
physical  processes  that  contribute  to  the  signal.  The  first 
term,  p(v",  represents  the  signal  that  is  related  to  the 
probe-induced  spatial  modulation  of  the  total  atomic 
population.  Such  a  contribution  has  not  been  considered 
before  in  similar  contexts.  As  we  will  see  below,  it  is  a 
signal  due  to  the  recoil-induced  resonances  involving 
both  the  probe  and  the  pump  fields.  The  term  p*"  deter¬ 
mines  the  nonvanishing  component  of  the  ground-state 
orientation,  while  p|J 1  and  p',°’_ ,  determine  different  com¬ 
ponents  of  the  atomic  ground-state  alignment.  In  the 
limit  that  |A|  »T,  they  represent  the  signal  due  to  Ra¬ 
man  processes  between  the  atomic  ground-state  sublevels. 
Both  contributions  are  analyzed  in  the  next  section. 

III.  CALCULATION 

OF  THE  PROBE  ADSORPTION  SIGNAL 

In  this  section,  we  derive  in  detail  the  probe  absorption 
coefficient  given  by  Eq.  (14).  Our  calculation  is  based  on 
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the  generalized  optical  Bloch  equations  described  above.  terms  proportional  to  x'  which  do  not  contribute  to  p ge 
In  the  low-field-intensity  and  low-energy  approximations,  to  linear  order  in  probe-field  strength),  the  equations  for 
and  keeping  terms  to  first  order  in  only  (and  dropping  the  ground-state  density-matrix  elements  are  given  by 


P*  1.*  i (/»./>')=-  r«+iP1JL 


"WFTA)^1-' -'{p  +  2*k'p,)-^^P* -^tp.p'  +  2*k)  < 

+  rf  dqN  (q)pg  Ug  ^p  +  fik +fq,p' +  fik +fiq) 

+  [  dq  N0(q]pg{)gl)(p+fk+fq,p'  +  fik+fq)  & 
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/ 

j  ntis  cra&i 

i  UT1C  TAB 

0 

□ 

i  Uii'inirOr  'coa 

□ 

1  By 

Distribution  / 
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Codes 

ciO/Of 

p,  '  -i  -«-c  lai 

+  —  f  dq  N  +  {q)[pg  _  |(p  — M  +Hq,p'—Hk  +  fiq)+pgi  giip  +  fik  +  fiq,p'+fi k  +  fiq)  t 

+pg  ,  gl(p -fik  +  fiq,p'  +  fk  +  fq  )+pgi ,g  \ip  +  fk  +  fiq,p'-fk  +  fiq  )]  1 
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+  ~~VTT->:  [dq  N^q)p  0  n(p +fk +*q,p' -fik  +  /tq) 
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rg+iP2 Pg'.g\iP’P’)+rA°Li{P’P')-irPgw{P’P') 

+  rf  dq  N+(q)pgt  gt{p-fik  +fiq,p' -fik  +fiq)+  ~  f  dq  N0(q)pg0g0(p  —  fik+fiq,p‘  —fik  +  fiq) 

pi 

+  —  J  dqN_(q)[pg_lg_l(p-fik+fiq,p'  —  -Fik+fiq)+pglgl{p  +  jik+fiq,p'  +  fik+fiq) 
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In  the  above  equations.,  yg  is  an  effective  ground-state  de¬ 
cay  rate,  which  could  be  determined,  for  example,  by 
transit-time  effects,  and  are  the  values  of 

Pgm.gm  (P’P')  t0  zeroth  order  in  probe-field  strength.  In 
our  calculation,  we  assume  that  the  pump  fields  have 
created  an  equilibrium  state  before  the  addition  of  the 
probe  field.  Therefore  ^g„,sm’(p,p')  are  the  equilibrium 
values  of  the  atomic  density-matrix  elements  in  the  pres¬ 
ence  of  the  two  pump  fields  only.  If  the  atomic  density 
matrix  is  diagonal  in  momentum  space  before  the  pump 
fields  arc  turned  on,  then  it  is  easy  to  deduce  from  the 
above  equations,  when  =  0,  that  the  only  non  vanishing 
density  matrix  elements  are  given  by  [3] 

pj,1 m.gm(P’P)=Pgm.gm(P >>  m  =  -l,0,l,  (24) 

and 

ft?U.g  .  ,(p±M,p  +*k  )=ff^Utn (p)  ■  (25) 

Note  that  the  atomic  momentum  distribution  function  in 
the  absence  of  the  probe  is  defined  by 

pt°W>  =  i  F?^.gJp,p)-  (26) 

m  =  I 


A.  Recoil-induced  signal 

To  calculate  the  recoil-induced  signal,  one  has  to  ob¬ 
tain  the  solution  for  ps(p,p  +2ftk  )  [see  Eq.  (14)].  In  or¬ 
der  to  analyze  the  effects  of  the  pump  fields  on  the  signal, 
it  is  convenient  to  use  the  momentum-family  notation  in¬ 
troduced  in  Ref.  [3],  For  example,  the  ground  states  that 
belong  to  the  same  family  J(p)  are  \gm,p+mfik), 
m  =  —  1,0,1.  The  stimulated  processes  involving  the 
pump  fields  do  not  change  the  momentum  of  a  given  fam¬ 
ily.  However,  the  addition  of  the  probe  field  can  induce 
stimulated  transitions  between  different  momentum  fami¬ 
lies,  a  process  that  leads  to  the  recoil-induced  resonances. 
Later  in  this  paper,  we  will  discuss  another  situation 
where  the  stimulated  transitions  involving  the  probe  and 
the  pump  fields  do  not  change  the  momentum  of  a  family. 

In  the  momentum-family  notation,  we  define  the 
ground-state  density-matrix  elements  that  are  related  to 
the  recoil-induced  signal  as 

pm(p)=pgmgm(p  +  rnitk,p  +  m*ik+2fik)  , 

m  =  -1,0,1  (30) 

and 


In  this  paper,  we  assume  that  the  width  of  atomic 
momentum  distribution,  which  is  of  order  p0  —  mu , 
satisfies  the  following  condition: 


(Ok 


«  ku  = 


kpo 


m 


«r . 


(21) 


The  above  inequality  implies  that  p0  »fik. 

Finally,  N .  n(q)  in  Eqs.  ( 1 9 )— ( 23 )  are  the  probability 
densities  for  emission  of  a  spontaneous  photon  having  a 
certain  polarization  (cr  or  tt)  and  a  z  component  of 
momentum  equal  to  Hq.  They  are  given  by  [7] 


A’,  (q 


1  + 


7 

11 

k* 


(28) 


In  the  present  calculation,  the  results  are  insensitive  to 
the  detail  of  the  functions  A \(q)  (e  =  ±0)  as  long  as  they 
satisfy  the  relation 

f  q.\tiq)dq  =0  ,  (29) 


which  they  do.  In  the  following,  we  use  a  single  normal¬ 
ized  probability  function  N(q),  which  also  satisfies  Eq. 
(29),  in  describing  the  momentum  distribution  of  a  spon¬ 
taneous  photon. 

_ I 


P  i.  -  i(P)=Pgt  i.s  *  i (p±*ik,p  +  fik  +2 ftk  )  .  (31) 

Furthermore,  the  ground-state  population,  and  orienta¬ 
tion  and  alignment  components  are  similarly  written  as 

Ps</>>=  2  pjpi  - 

m  -  -  1 

Por(P>=-^[p|(P)_P  1(P>]  -  ,32) 

Pai(P>=7T[Pi<P'+P  -  |(P)~2p0</>)]  • 
v  o 

Since  the  only  nonvanishing  density-matrix  elements  to 
zeroth  order  are  given  by  Eqs.  (24)  and  (25),  the  terms  in 
the  above  definitions  are  all  first  order  in  the  probe-field 
strength.  We  have  omitted  the  superscript  (1)  for  con¬ 
venience.  The  recoil-induced  signal  is  now  proportional 
to  the  integral  of 

f  p.ip'idp  —  J ps'  ’( p,p  +2fik  )dp  , 

as  is  evident  from  Eq.  (14).  We  now  proceed  to  obtain 
the  equation  for  the  evolution  of  ps(p). 

First,  p .  ,  .  ,(p)  are  obtained  from  Eqs.  (22)  and  (231  as 


P  u(Pl  = 
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5r' 


L_ 

y  —  i  &  1  ‘  y  +  i  A 


P  i<P> 


+  3  f  dq  N(q)[p  ,(p  +ftk  +ftq  )+pPp  -fik  +fiq  )  +  4po<P  *] 

"  rnT*ffi#’  u"  ~M-p  *m  ’■  7r(^7iTpJ>'"  m^*«M**»  ' 

+  7^^7777  f  ^W^pT  l.lti<p+fy,p  +  2ftq+fiq)  ■ 

5T  (y  +  A  ) 


(33) 


and 
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5P' 


y  —  i  A 


y+iA 


Pi<P> 


+  { f  dq  N{q)[p_l(p+fik+-fiq)+pl(p—  fik+fiq)  +  ±p0{p+fiq)] 


XX’*eib' 


5  r 


g_,(p+M,p+^fc)  Pgi  ,xi(p+^fc,p+^) 

y  —  /  A  y  +  /  A 


~<0) 


p",;g_j(p  +  3flk,p  +yi/c )  +  pgl  g _{(p+ftk,p-fik) 


y  —  i  A 


y +  2  A 


+  ■ 


E^!i 


•  _i6f 


[P^ilgifp  +  fiq ,p  +  ^  )  +^g°i i,g  - i(p  +%,P  +  % ) 


5T'(y2  + A2) 


*-!.*- 


+tfgi,g  -\(p  +  2fik  +fiq,p+fiq)  +  jftgo,g0(p  +tiq,p+fu} )  ]  , 


3*<0) 


(3' 


where  we  have  neglected  yg,  8,  and  the  kinetic-energy  term  (p2—p,2)/2mfi  in  comparison  with  P  [see  inequality  (27)]. 
Substituting  Eqs.  (33)  and  (34)  into  Eqs.  ( 19)  — (21 ),  and,  assuming  the  secular  limit 

r  / 1  a  | « i  ,  <3 

to  lowest  order  in  T /|  A|,  one  obtains  the  equation  for  ps(p)  as 


i>s(p)=  r’ps(p)  +  r/ dt 1  #(?){  ir^P^P  +flk  +ful  )+Pi<P -*k  +*<1  >1 


+  rf  dq  N(q) 


22^2 
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V2 

180 

Vi 

540 


+  ^[p^p  +2fik  +fiq)+  ps{p  -2iik  +fiq )]  +  ±ps(p  +fiq)\ 
[pjp  -fik+fiq) ~por(p  +M  +fiq)) 

tPor<P  -2  fik+fiq)  ~por(p  +2  tk+fiq)] 


-  rpjp )+~~rfdqN(q)\ pjp  +  2 M  +  *q )  +pjp  - 2 fik  +  fiq  ) 

+  44[pal(p  +ftk  +fiq  )+pa!(p  —  fik  +  fiq)]  —  66piX(p  +fiq)  j 


.  2  kp 

—  1 — r-  —  1 


i4 (Ol 
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ps(p)~i4akpot(p) 


i.« - 1 {P  +*k'P  +fik)-Xhtr^g - tP-*k'P -*k ) 


+  Fy2  +  ^r  Id<l"{(1)ti0)-lg-'{P+*<l'P+*(l) 


Xre'bl  -AO) 


XX'*e'b'  ^o, 


IfcrirfZJP- +2«k,P  +  2« 


+  ^2+^2)-  f  d<*  m^l^P  +*k  +*<l,P+*k  +fiq ) 

-  + 3*  •" + 5“ +«./>+«» 

+  f  dq  N(q)f?g°i\gl(p  +  2ftk+fiq,p+2*k+j>q) 

_  xarv-tfl UI(p -Hk,p +fik )- ,.g,(p +«.P +3« ) 
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J  dq  N(q)pglX  gX(p+fiq,p  +  2fik+fiq)  . 
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Note  that,  neglecting  the  photon  recoil  momentum  fik 
and  fiq ,  one  has 
ps(p)  =  0  . 

In  other  words,  when  the  photon  recoil  momentum  is 
neglected,  the  population  modulation  term  ps(p)  remains 
zero  despite  the  addition  of  the  probe  field.  This  is  in 
agreement  with  the  semiclassical  theory  where  the  atomic 
center-of-mass  motion  is  described  by  a  classical  variable 
(the  velocity),  and  the  total  atomic  population  is  a  con¬ 
stant  [5],  In  this  case  the  population  will  remain  diagonal 
in  momentum  space  before  and  after  the  addition  of  the 
probe  field  [see  Eq.  (24)].  It  is  clear  here  that  the  signal 
related  to  the  probe-induced  modification  of  p5(p)  is  a 
direct  result  of  the  photon  recoil  momentum,  which  ex¬ 
plains  why  we  call  this  part  of  contribution  to  the  probe 
absorption  the  recoil-induced  signal. 

From  Eqs.  (19 )— ( 21),  we  then  solve  for  the  quasi¬ 
steady-state  solutions  for  por(p)  and  pal(p)  in  terms  of 
ps(p),  and  substitute  the  solutions  into  Eq.  (36).  In  the 
secular  limit  (35),  por(p),  whose  value  is  proportional  to 

_ 1 


(  T /A )2{(0k  /ku  )ps(p),  can  be  neglected,  while  pa|(p)  is  ap¬ 
proximately  given  by 

P°'{P)=^rP*{P)  '  (37) 

Corrections  to  Eq.  (37)  are  smaller  by  a  factor  ojk  /  ku). 
We  have  dropped  some  terms  in  the  solution  for  pa,(p) 
that  are  explicity  proportional  to  ¥'*■  Since  their  contri¬ 
butions  to  the  recoil-induced  signal,  upon  substitution  of 
pal(p)  into  Eq.  (36),  is  of  order  co2k  /(ku)1  [8],  and  since  we 
are  interested  only  in  the  signal  that  is  linear  in  uk  /( ku  > 
[2],  such  x  *  terms  can  be  neglected  when  replacing  pjp) 
in  Eq.  (36)  with  Eq.  (37).  After  such  procedures,  we  ar¬ 
rive  at  an  equation  involving  p4(p)  only.  In  this  equation, 
there  are  two  types  of  terms.  First,  there  are  terms  which 
are  explicitly  proportional  to  \  *  ■  They  are  identical  to 
those  appearing  in  Eq.  (36),  and  provide  the  source  terms 
for  the  recoil-induced  signals.  Second,  there  are  terms 
which  are  proportional  to  T'.  They  describe  the  optical 
pumping  effects  due  to  the  presence  of  the  pump  fields, 
and  are  given  by 


[p, <P > ]pump  =  -  f§5  r'P, <  A ' +  rf, dq  N(q)[p, (p  +2*k+*q)+px(p-  2M  +*q)] 

+  —jT' J" dq  A r(q)[ps(p  +ftk  +fiq  >+p,(p  —  ftk+fiq)]  +  T'  f  dq  N(q)ps(p  +fiq)  .  (38 1 


If  one  assumes  that  N(q)  —  b(q)  and  expands  the  in¬ 
tegrands  around  p,  it  can  be  simplified  as 

[Av<P>]pump=T?r<^  )2^TPS(/,)  ’  (39) 

dp- 

where  the  overdot  symbol  in  Eq.  (39)  denotes  time  deriva¬ 
tive  only.  Equation  (39)  describes  a  momentum  diffusion 
process  induced  by  the  pump  fields  in  the  absence  of  the 
cooling  forces  (since  we  have  neglected  the  Doppler  shift 
in  comparison  with  T').  The  diffusion  coefficient, 
A V(fik  )2,  is  in  agreement  with  the  result  of  Ref.  [3]  in 
the  limit  |A|»r  and  neglecting  the  contribution  of 
spontaneous  emissions. 

In  this  paper,  we  are  interested  only  in  the  qualitative 
effects  of  the  optical  pumping  on  the  recoil-induced  sig- 
_ 1 


I  -  " 

nal.  For  the  sake  of  simplicity,  and  without  loss  of  the 
basic  features  involved,  we  replace  the  optical  pumping 
terms  in  Eq.  (38)  by  the  following  simplified  model: 

[  As  <  A  ]Pump=_rp[2P«<P>_P.v  </>+*»*  ]-pSp-*k  I]  • 

(40) 

where  T p(fik )2  is  the  atomic  momentum  diffusion 
coefficient  due  to  the  pump  fields,  and  rp  is  of  order  T’. 
One  can  verify  that,  by  expanding  thep,(p±M  )  terms  in 
the  above  equation  around  p,  Eq.  (40)  leads  to  an  equa¬ 
tion  for  momentum  diffusion  that  has  the  same  structure 
as  Eq.  (39).  The  complete  equation  for  p5(p)  can  now  be 
written  as 


Av<A  =  -rP[2p,(p'-pv(p+M  )~ps(p—fik  )]- 


P«<A 


u  xlp+*k.p+*k)~  I.*  ,<P-**.P-#*) 


+  fd^N{^]pKh  U  ttp+fiq,p  +  hq) 


pfA 

„  ibt 


yy'* 


2(  Y 


+  f  dq  N(q) fi'^p +*k  +*q,p+*k  +fiq) 

_  XX"^'  ,'<) 
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'•-'Of  yy»*  ,)0/ 

-_7a]  +  m -P  +  ™k)~  bX{LiA' /"Up  +  ^k 
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+  f dq  N^(p  +  2Hk  +*q,p+2*k  +*q  ) 

6(y‘  +  A*) 


y  y  ^  ®  y  y  ^  ^ 

>  *  „  itu 


+  fdqmq)ft°lx  ^p+*q,p+  2M  +  liq  )  . 

6(y  +  A*' 


(41) 


It  is  shown  in  Appendix  B  that  by  imposing  the  following 
condition: 


£'  (t>l 
ku  ( ku )2 


(42) 


one  can  neglect  the  effects  of  the  optical  pumping  terms 
in  Eq.  (41).  Since  the  momentum  diffusion  coefficient  is 
of  order  T'rfr  /k1,  and  the  atomic  momentum  distribution 
width  is  of  order  mu,  the  above  condition  implies  that, 
during  the  time  of  order  \/ku  in  which  the  ensemble  of 
p,°\p,p  +2 ftk  )  decays  as  a  result  of  Doppler  dephasing, 
the  amount  of  momentum  diffusion  due  to  optical  pump¬ 
ing  is  small  compared  with  the  initial  momentum  width 
mu.  Assuming  that  condition  (42)  is  satisfied,  the  terms 
on  the  right-hand  side  of  Eq.  (41)  that  are  proportional  to 
T,  can  be  dropped.  Fn.ally,  we  assume  the  zeroth-order 
terms  in  Eq.  (41)  are  given  by  [3]. 

$'u-tp)=p'?\mP**k) , 

/W<P>=Po0,^(/>)  ,  (43) 

pT  U  •  i<P)=Pm.  •  1  W(p)  , 


where  W{p)  is  a  normalized  distribution  function  of  p 
which  has  a  width  of  order  p0~mu  that  satisfies  inequali¬ 
ty  (27). 

The  population  modulation  term  ps(p)  can  now  be 
solved  from  Eq.  (41)  without  the  optical  pumping  terms, 
and  integrated  over  p  to  give  the  recoil-induced  signal  to 
lowest  order  in  iok  /(ku)  and  T/|A|  as 


___  JO  vVy'*IyI2 
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(44) 


The  function  I (x)  is  defined  as 


I(x)~ 


Vi 


f 


IV(pnt  )pndt 


(45) 


and  /'(x)  is  the  first-order  derivative  of  /(x)  with  respect 
to  x.  The  values  of  p'£'  can  be  deduced  from  the  results 
of  semiclassical  calculations  [3].  In  the  limit  that 
ku  «  r\  they  are  given  by 


pT- 


-p{,)\- 


n 
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„ini_ 
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p  1.1 


I - — 

As  shown  in  Ref.  [2],  the  signal  Im[p£,  ]  as  a  function  of  5 
has  a  dispersion-like  line  shape  around  5  —  0  for  a  Gauss¬ 
ian  distribution  of  H/(/>),  and  it  has  the  same  sign  as  we 
have  predicted  qualitatively  in  the  Introduction. 


B.  Raman  signal 


We  now  consider  !he  contribution  to  the  probe  absorp¬ 
tion  signal  due  to  the  atomic  ground-state  orientation  and 
alignment,  which  can  be  interpreted  as  a  Raman  signal  in 
the  limit  that  |A|/T  »1.  Since  both  the  ground-state 
energy  splitting  !S£i  -  iY!2lA|/(y2  +  A2)  and  the  effective 
decay  rate  T'  of  each  sublevel  are  much  larger  than  iok , 
one  can  neglect  the  recoil  shift  terms  in  considering  the 
Raman  signal.  In  this  case,  the  calculation  for  the  Ra¬ 
man  signal  can  be  carried  out  semiclassically  in  a  sense 
that  the  atomic  center-of-mass  momentum  and  position 
are  treated  as  classical  variables.  The  semiclassical  equa¬ 
tions  for  the  density  matrix  can  be  obtained  directly  from 
Eqs.  ( 19)-(23)  by  the  transformation  [2] 


{zJ)-^ffdPdP'exP 


ft 


Xp^^p.p')  ,  (47) 

neglecting  terms  proportional  to  the  recoil  energy  fuok. 
In  Eq.  (47),  j  is  the  classical  atomic  center-of-mass  posi¬ 
tion,  which  can  depend  on  time.  Defining  similar  quanti¬ 
ties  for  the  atomic  ground-state  population,  orientation 
and  alignment  as 

i 

P\  2  Pxm.gm(Z,t)  , 

m  l 

Pot  IPkI.kI^2’^  Px  l.Jf  1^2,^1  ’ 

(48) 

Pi||=  ^[Ps-'.Sl(z',)  +  Pv  U  1(2’'> 


~2Px ruiM'M  • 


P  i.  -  i=Pri.«  •  i*2-'*  ’ 


(46) 


one  then  obtains  the  following  equations. 
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As  is  evident  from  Eq.  (49),  the  total  population  in  the  where  r  is  the  z  component  of  the  atomic  velocity.  One 


semiclassical  picture  is  a  conserved  quantity,  i.e., 
Ps=  1  • 


can  verify  that  to  zeroth  order  in  kv /Y\  Eqs.  (54)  and 
(56)  lead  to  the  results  of  Eq.  (46)  in  Sec.  Ill  A. 

Upon  substitution  of  Eq.  (56)  into  Eqs.  ( 50 )  —  ( 53 ),  one 


We  now'  make  the  following  e-  p  ...sions  in  powers  of  the  can  then  solve  for  the  ground-state  elements  to  first  order 
probe-field  strength  for  the  steady-state  solutions  of  Eqs.  ’n  probe-field  strength.  In  the  limit  that  ku  /T  "0,  the 
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first-order  terms  are  given  by 

p.i»=_ivJ. _  rv _ 

/,,r  204  r'/6  +  /5  +  (r'/3):( A/D:/(5r'/6  +  /6) 

x  3_ _ 5r*/3 _ iA  1 

5r'/6  +  /8  T  y  —  i  A 


where  the  superscript  (0)  or  (1)  again  indicates  that  the  j'  +  t(A  +  8) 

related  terms  are  to  zeroth  or  to  first  order  in  probe-field  p|  j'  =  o  , 
strength.  Equations  (  50)- ( 53 )  can  then  be  solved  recur-  ^ 

sively.  First,  the  zeroth-order  terms  are  solved  by  letting  p\"  — — - p\ 

\’  =  0  in  Eqs.  ( 50 )— ( 53 ).  One  finds  3v  3  ^  5T /6  +  /ft 
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The  Raman  signal  [see  Eq.  (14)]  is  a  function  of 
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Upon  substitution  of  the  results  (57)  into  Eq.  (58).  one  ob¬ 
tains  the  Raman  signal  to  zeroth  order  in  ku  /V  as 


47 


RECOIL-INDUCED  RESONANCES  IN  PUMP-1  ROBE  .  . 


4130 


-4  Rm » _ 

Pc.'  - 


25 

1224 


<v"  v  : 

r'/is  /A 
5r/6+/8  r 

y  +  i  A 

rV6  +  i8+,r^)2«A/[>2 
5T  /6  +  /8 

3 


In  the  limit  A  >  >  T,  the  Raman  signal,  proportional  to 
Imliv'"  ]•  exhibits  resonanees  at 

b  •  •  ^  y  ■  (60) 

Both  resonances  are  Lorentzian  type  with  opposite  signs, 
and  the  full  width  at  half  maximum  (FWHM)  of  each  res¬ 
onance  is  given  by  F'.  Ft  negative  A,  the  probe  gain 
and  probe  absorption  occur  when  8  equals  r'A/(3D  and 
r'A/(3D,  respectively.  The  ratio  between  the  magni¬ 
tudes  of  the  gain  and  the  absorption  peaks  is 

The  above  result  can  be  understood  if  one  chooses  the 
quantization  axis  along  the  local  direction  of  the  total 
pump-field  polarization.  As  shown  in  Ref.  [3],  the  com¬ 
bination  of  the  counter-propagating  a  *  —  cr~  pump  field 
results  in  a  linearly  polarized  field  whose  strength  is  space 
independent  and  whose  polarization  vector  rotates  along 
the  laser  propagation  axis.  Relative  to  the  local  quantiza¬ 
tion  axis,  the  pump  field  is  jt  polarized  with  an  amplitude 
equal  to  v/26  . ,  while  the  probe  polarization  has  a  '  ,o  , 
and  tt  components.  Owing  to  different  Clebsch-Gordan 
coefficients  for  different  transitions  (see  Fig.  2),  to  zeroth 
order  in  probe-field  strength,  the  steady-state  population 
of  the  m  =0  state  is  greater  than  those  of  the  m  =±1 
states,  and  the  light  shifts  of  these  states  induced  by  the 
pump  field  also  differ.  One  can  verify  easily  that  the  en¬ 
ergy  splitting  betw  “n  the  m  =0  and  ±1  states  60  . ,  is 
given  by  Eq.  (60).  For  |A|»r  and  negative  A,  when 
5<0,  the  dominant  process  is  the  transition  from  m  =0 
to  .■*  1  states  via  absorption  of  a  pump  photon  and  emis¬ 
sion  of  a  probe  photon,  resulting  in  probe  gain;  when 
8>0,  the  same  transition  involves  absorption  of  a  probe 
photon  and  emission  of  a  pump  photon,  resulting  in 
probe  loss.  As  shown  in  Fig.  2(a),  these_  two  pro- 
cesses  _  have_  relative  strength  [( 1  /  v^ 6IV7 2/3]2/ 

[( 1  /v/2H  l/t/2))2  — £. 


C.  Overall  signal 


The  overall  signal  is  proportional  to 

Im(p,,.)  =  Im[p/;’]  +  Im[/^rT^  ,  (61) 


which  is  a  superposition  of  the  recoil  and  Raman  signals. 
When  IAI  »y,  both  signals  are  proportional  to  1/|A|2. 
The  dimensionless  parameter  that  determines  the  magni¬ 
tude  of  the  recoil  signal  to  that  of  the  Raman  signal  is 
given  by 


Mk  r 
tiku)1 


(62) 


where  £  <  1  is  a  parameter  that  characterizes  the  popula- 


5rV3  /A 

5rv6+/8  r  t  2  _ 

y  —  i  A  y  +t(  A  +  8 ) 


(59) 


+  ■ 


1 


1 


5r'/6  +  /6  y-iA 


I  j 

tion  difference  between  the  initial  and  final  states  of  the 
Raman  transition. 

Figures  3(a)  and  3(b)  show  the  probe  absorption  signal 
cs  a  function  of  8  for  two  different  sets  of  parameters, 
with  different  values  of  A/y.  The  recoil-induced  signal  is 
centered  at  8  =  0,  having  width  of  order  2ku.  The  Raman 
signal  has  a  width  equal  to  T'.  The  distance  between  the 
two  Raman  peaks  corresponds  to  the  differential  light 
shifts  between  different  magnetic  substates.  As  one  can  i 

see  from  Fig.  3,  when  !  A!  /}'  increases,  the  splitting  of  the 
two  Raman  peaks  also  increases,  while  the  recoil-induced 
signal  remains  a  dispersionlike  curve. 

For  ideal  sub-Doppler  cooling,  one  can  utilize  a  rela¬ 
tion 


-0.20  -0  10  0.0  0  10  0  20 
6/y 


FIG.  3.  Probe  absorption  coefficient.  The  dashed  line  shows 
the  contribution  of  the  Raman  resonances,  while  the  dot-dashed 
line  shows  that  of  the  recoil-induced  resonances.  The  solid  line 
shows  the  overall  signal.  For  both  figures,  ti>t  /y  =  0.0008,  and 
y„/y  =0.0001.  (a)  A/y=-5.0,  y/y=I.O,  ku/y  =  0.02;  (b) 
A/y  =  -  10.0,  \/y  =  3.0,  ku/y  =0.04. 
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tok  |A 

( ku ):  kHT  |^|: 


(63) 


to  rewrite  >/  as 


For  values  of  iAi  /}'  that  are  of  order  1  /£,  the  magnitude 
of  the  recoil-induced  probe  absorption  signal  is  compara¬ 
ble  to  that  of  the  Raman  signal.  Equation  (63)  also  indi¬ 
cates  that  the  width  of  the  recoil-induced  signal,  which  is 
of  order  ku,  increases  linearly  with  !  y|  for  a  fixed  value  of 
!  A  i .  This  is  in  contrast  with  the  Raman  signal,  whose 
width  [~T']  increases  quadratically  with  |y|.  Finally, 
when  the  cooling  condition  is  not  optimal  and  ku  is  of  or¬ 
der  F',  the  Raman  signal  will  be  a  convolution  of  a 
Lorentzian  function  of  width  T'  with  a  Gaussian  function 
of  width  1 .66(  2A«  ). 


IV.  DISCUSSION 

We  have  presented  the  calculations  of  probe  absorption 
coefficients  in  the  JK  =  1  — <-Je  =2  systems  in  the  presence 
of  polarization  gradient  pump  fields.  Apart  from  the  Ra¬ 
man  signals  which  result  from  the  nonconservation  of  the 
atomic  properties  such  as  the  ground-state  orientation  or 
alignment,  there  is  a  recoil-induced  signal  which  is  relat¬ 
ed  to  a  probe-induced  spatial  modulation  of  the  atomic 
total  population.  We  find  that  optical  pumping  does  not 
broaden  the  recoil-induced  resonances.  This  can  be  un¬ 
derstood  as  follows.  In  the  picture  of  Raman-type  pro¬ 
cesses  between  center-of-mass  momentum  states,  the 
transition  width  is  determined  by  the  effective  lifetime  of 
each  momentum  subclass.  When  condition  (42)  is 
satisfied  [(  V /ku)Wo\  /k  2«-)«l],  the  population  in  a 
given  momentum  interval  A p  does  not  decay  during  the 
time  interval  1  /(ku).  As  a  result  the  optical  pumping 
does  not  broaden  the  linewidth  of  the  recoil-induced  sig¬ 
nal. 

The  recoil-induced  signal  calculated  in  Sec.  Ill  is  relat¬ 
ed  to  the  stimulated  transitions  between  different  momen¬ 
tum  families  involving  the  pump  and  the  probe  fields,  ow¬ 
ing  to  this  particular  probe  configuration.  On  the  other 
hand,  if  the  probe  is  a  a*  field  copropagating  with  the 
o  pump,  the  stimulated  transitions  involving  the  probe 
field  are  between  states  that  belong  to  the  same  family. 
As  a  result,  there  is  no  recoil-induced  signal  in  this  case. 
However,  it  has  been  shown  recently  by  Lounis  el  al.  [9] 
that  there  is  a  Rayleigh-type  resonance  in  this  situation 
even  without  the  inclusion  of  recoil  effects,  owing  to  some 
nonvanishing  spatially  averaged  atomic  velocity  that  os¬ 
cillates  at  a  frequency  6  under  the  influence  of  the  copro¬ 
pagating  pump  and  probe  fields  having  the  same  polariza¬ 
tion.  This  average  velocity  leads  to  a  so-called  motion- 
induced  atomic  orientation  [3],  which  scatters  photons 
from  the  pump  field  into  the  probe  field,  resulting  in  the 
probe  gain  or  loss.  Such  a  Rayleigh  resonance  will  not  be 
broadened  by  the  atomic  Doppler  width  since  it  only  in¬ 
volves  transitions  between  atomic  states  of  the  same 
momentum.  It  is  clear  from  the  above  arguments  that 
there  is  no  Rayleigh  resonance  for  the  field  configuration 


considered  here,  since  the  spatially  averaged  atomic  ve¬ 
locity,  to  first  order  in  the  probe-field  strength,  is  zero. 
As  a  result  there  is  no  motion-induced  atomic  orientation 
in  this  case.  Instead,  there  are  recoil-induced  resonances, 
i.e.,  Raman  processes  between  atomic  states  of  different 
kinetic  energies,  for  this  field  geometry.  The  width  of  the 
recoil-induced  signal  is  determined  by  the  Doppler  wadth 
2 ku  and  serves  as  a  probe  of  the  atomic  velocity  distribu¬ 
tion. 

The  theoretical  calculations  presented  here  cannot  be 
compared  directly  with  the  experimental  results  on 
pump-probe  spectroscopy  of  atoms  confined  in  a 
magneto-optical  trap  [4].  One  reason  is  that  the 
actual  atomic  level  schemes  are  more  complicated 
(F=4— *F  =  5  for  Cs  atoms)  than  that  considered  in  this 
work.  Also  the  possible  inclusion  of  inhomogeneous 
magnetic  fields  in  those  earlier  works  can  further  compli¬ 
cate  the  problem.  In  the  results  of  some  more  recent  ex¬ 
periments  in  one-dimensional  (ID)  optical  molasses  pro¬ 
duced  by  a  pair  of  circularly  polarized  fields  [9,10],  there 
have  been  some  signatures  of  the  recoil-induced  reso¬ 
nances  considered  in  this  and  a  previous  paper  [2].  How¬ 
ever,  a  direct,  unambiguous  observation  of  such  reso¬ 
nance  phenomena  is  yet  to  be  made.  When  selecting  the 
field  configuration  as  the  one  considered  here,  it  may  be 
useful  to  choose  parameters  such  that  |  A  i  / T  » 1 ,  there¬ 
by  the  normal  Raman  signals  can  be  separated  from  the 
recoil-induced  signal  centered  at  6  =  0. 

Finally,  we  discuss  the  relation  between  the  results  cal¬ 
culated  here,  and  a  number  of  experiments  on  cold  atoms 
using  linilin  cooling  field  configuration  [11,12],  in  which 
the  linear  polarizations  of  the  counterpropagating  cool¬ 
ing  fields  are  orthogonal.  As  is  well  known,  atoms  in 
such  field  configurations  experience  an  effective,  spatially 
modulated  potential  originating  from  the  atomic 
ground-state  energy  shifts  in  the  fields.  If  the  average 
atomic  kinetic  energy  is  cooled  below  such  potential 
depth,  a  significant  portion  of  the  atomic  population  be¬ 
comes  spatially  localized  [13].  Such  localized  atoms  have 
a  bandlike  energy  structure  [14],  which  can  be  probed  by 
various  spectroscopic  means.  In  situations  like  this, 
where  spatial  localization  of  atoms  is  important,  the 
recoil-induced  resonances,  which  are  basically  a 
phenomenon  associated  with  the  atoms  in  the  energy 
continuum  states,  may  still  exist.  Due  to  the  small  popu¬ 
lation  of  untrapped  atoms  in  most  situations,  and  also 
possibly  due  to  the  large  amplitude  of  the  signal  provided 
by  the  trapped  atoms  [15],  the  recoil-induced  signal  may 
be  too  small  to  be  observable  under  most  experimental 
conditions.  There  might  be  some  cases,  for  example,  with 
low  cooling  field  intensities,  where  the  localization  effects 
arc  less  significant,  and  the  recoil-related  effects  are  more 
pronounced.  Whether  such  situations  exist  requires  fur¬ 
ther  research. 
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APPENDIX  A 

In  this  Appendix,  we  give  the  derivation  of  the  probe 
absorption  coefficients  in  terms  of  the  atomic  ground- 
state  density-matrix  elements.  A  convenient  basis  for  the 
derivation  is  the  irreducible  tensor  basis  [5,16],  in  which 
density-matrix  elements  are  defined  as 

p$(F,,F2)=  2  (-\)Fl~m2(Fl,mi;F2,-rn1\K,Q) 

m  j ,  m 

'  ^1) 

where  -f  m  are  density-matrix  elements  written  in 
the  magnetic-state  basis,  and  (Fl,mx;F2—  m2\K,Q)  is  a 
Clebsch-Gordan  coefficient.  The  polarization  of  the 
medium  can  be  written  as  [16] 

P  =  -jT(-l)H~c  +  1er(,g  2  eqp\(G,H)+ c.c.  ,  (A2) 

1  $  9  =  ±1.0 

where  G,H  are  the  angular  momenta  of  the  ground  and 
excited  states,  respectively,  reg  is  a  reduced  matrix  ele¬ 
ment  of  the  dipole  moment,  and 

e,i=  + -4|-(x±/y),  e0=z  .  (A3) 

The  absorption  coefficient  for  a  probe  field  having  (com¬ 
plex)  polarization  e',  frequency  O',  and  propagation  vec¬ 
tor  k',  is  proportional  to  [17] 

e,‘-P(k')==2(”1)?<jP-9(k') 

9 

er 

=  -^(-l),/-G  +  12<pJ(G,tf;k')  ,  <A4) 

3  9 

where  P(  k' )  and  pq(G,H; k' )  denote  the  components  of  P 
and  pq(G,H)  that  vary  as  exp(  —  /k'-R  +  'ft'f ),  and  the 
e^’s  are  given  by 

*'•  j  =  +  =  <-  .  <A5) 

where  e'x,  e'v,  and  e'z  are  the  Cartesian  components  of  the 
probe-field  polarization  vector. 

In  the  case  of  Sec.  Ill,  the  probe  polarization  is  given 
by 

C,  =  -&,.!•  (A6) 

Substituting  Eq.  (A6)  into  Eq.  (A4),  one  obtains 

pge*p\lG,H;k)  ,  (A7) 

where 


Pgc  Pg  -  l.t  -  2^  k  )  “b  ^PgO,e  - k) 


+  ^PgyJk)  e  iXl,  <A9) 

where  pgm,em‘lk)  are  related  to  density-matrix  elements 
Pgm.em'^P'P'  1  defined  in  Eq.  (9)  through 

Pgm.em-lk)=  f  Pgm,,m(P<P  +*k  )dp  .  (A10) 

When  replacing  p with  ground-state  matrix  ele¬ 
ments,  and  neglecting  terms  that  correspond  to  linear 
probe  absorption,  one  obtains  the  following  expression 
for  pgc- 

Pg,  =  ~  +  /  Ip-  i < P,P  +  2 M  )  +  \p<}  \p,p  +  2 fik  ) 


+  +2/tk  )]dp 


iX\e~ibt 


7 S  ftu-\[P'P)dP 


where  p}\  and  ptf  ‘  denote  pgl\  g.,  and  p^o  respective¬ 
ly.  By  using  the  identity 

P-  i  +  iPa+^Pi^liPs-  h(P\~P -t* 

+  T(p,+p..,-2p0),  (A12) 

where  ps  is  the  total  population  as  defined  in  Eq.  (16),  one 
obtains  Eq.  (14)  for  pgc. 


APPENDIX  B 


In  the  appendix,  we  analyze  the  effects  of  optical 
pumping  due  to  the  presence  of  the  pump  fields  on  the 
recoil-induced  signal,  and  give  the  derivation  of  condition 
(42).  The  optical  pumping  processes  cause  an  effective 
atomic  momentum  diffusion,  which  can  be  represented  by 
Eq.  (40).  Instead  of  using  an  effective  transit  decay  rate 
}'g,  we  solve  the  time-dependent  problem  for  the  recoil 
signal.  The  evolution  equation  of  p(p)  can  be  written 
from  Eq.  (41)  as 

ps(p)=  -2rpp,(p)  +  2rp/ <o(p'-+p)psip')dp' 

+  i  ^+40*  p{p)  +  e‘b'W0(p)  ,  (Bl) 

m 

where  H^tp)  represents  the  “source”  terms  in  Eq.  (41) 
(terms  explicitly  proportional  to  y'*).  The  kernel 
a>(p'  — *p )  is  given  by 

a)(p'~-*p)=  \[b(p‘ -p+ftk  )  +  8(p'—  p-fik  )]  .  (B2) 


Equation  (Bl)  can  be  solved  in  the  interaction  representa¬ 
tion  defined  by 


p,(p)=p(p)e 


i l  2kp  /m  4  )f 


(B3) 


p\[G,H\k)=p\(G,H,k)eihl  . 
In  the  magnetic-state  basis,  one  has 


Moreover,  by  noting  that  «(p'-*p)  depends  only  on  the 
difference  (p  —  p'),  one  can  make  a  Fourier  transforma¬ 
tion. 
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p(p)=  , 

and  obtain  the  solution  for/(£)  as  follows; 


(04)  centered  at  p—  0,  with  the  most  probable  momentum 
p0  =  mu.  Then  upon  substitution  of  Eqs.  (B5MB7)  into 
Eq.  (B8),  one  obtains 


MS)=exp[-2rpt  +  2rpf‘5)(£,t')dt'] 

X  J"Jexp[ 2rgt' ~2rp  f^fi>(t;,t")dt"]a)s{£;,t' )dt '  , 

(B5) 

where 

tu(£,/  )  =  yi  exp[  —ifiL,k  —  1 4ro*  t  ] 

+  exp[//ifA;+/4«jl/]|  ,  (B6) 

and 

.  (B7) 

The  recoil  related  probe  absorption  signal  depends  on 
fe  ,s'p(p),dp  [see  Eqs.  (14),  (15),  (30),  and  (32)].  From 
Eqs.  (B3)  and  (B4),  it’s  easy  to  show  that 


fe  ih'p(p)sdp=// 

-2k 
- 1 

m 

(B8) 


As  an  illustration,  let  us  assume  fV0(p)  to  be  a  Gaussian 


Jp[p)dp=  f'exp 


-2Vpt'  +  2Tp 


sin(4  io,t') 


4w, 


Xexp[  —  (kut')1 2 3 4 5 6 7 8  —  ibt'  +  i4a)kt']di‘  . 


(B9) 


It  is  obvious  that  the  values  of  t'  that  contribute  to  the 
above  integral  are  t'5  \/(ku).  All  the  efTects  of  the  opti¬ 
cal  pumping  are  included  in  the  first  exponential  term  on 
the  right-hand  side  of  Eq.  (B9).  In  order  for  the  pumping 
effects  to  be  negligible,  this  term  has  to  be  approximately 
unity.  Since  cok  /iku)  «\,  one  can  expand  the  sine  func¬ 
tion  and  reach  the  following  expression; 


exp 


4  rP  frk 2 
3  ku  m2u2 


(BIO) 


which  is  true  when  condition  (42)  is  satisfied.  The 
momentum  diffusion  coefficient  due  to  the  optical  pump¬ 
ing  effects  is  T pfi2k2.  Condition  (42)  implies  that,  during 
the  coherence  dephasing  time,  the  amount  of  momentum 
diffusion  due  to  optical  pumping  is  small  compared  with 
the  width  of  the  original  momentum  distribution. 
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